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Abstract 

We suggest kinetic models of dissipation for an ensemble of interacting ori- 
ented particles, for example, moving magnetized particles. This is achieved by 
introducing a double bracket dissipation in kinetic equations using an oriented 
Poisson bracket, and employing the moment method to derive continuum equa- 
tions for magnetization and density evolution. We show how our continuum 
equations generalize the Debye-Hiickel equations for attracting round particles, 
and Landau-Lifshitz-Gilbert equations for spin waves in magnetized media. We 
also show formation of singular solutions that are clumps of aligned particles 
(orientons) starting from random initial conditions. Finally, we extend our 
theory to the dissipative motion of self-interacting curves. 
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1 Introduction 

Many physical systems - such as plasmas, particle-laden fluids, self assembling molecules 
etc. - can be understood in terms of collective dynamics of a large number of particles 
moving in an ambient medium. Dissipation plays a crucial role for these systems, in 
forcing them towards an eventual equilibrium state in the absence of external forces. 
The dissipation in these heterogeneous systems is fundamentally different from inter- 
nal dissipation in homogeneous fluids, which is usually described by diffusion terms 
in momentum equations (like the Navier-Stokes equations). The goal of this paper is 
to demonstrate a unifying principle for deriving dissipation in systems consisting of 
particles which are embedded within an ambient fluid, or matrix. The only role of 
that matrix is to exert friction on individual particles at the microscopic level, and 
the matrix itself is assumed to have no dynamics of its own. This friction is applied 
both in configuration space and in momentum space, which leads to the necessity of 
considering kinetic equations. 

The physical picture we have in mind consists of microscopic magnetized particles 
moving under the influence of external forces and mutual attraction in a passive 
ambient continuum. Such systems have attracted much attention recently because of 
their technological significance. For example, magnetic nanoparticles in this setting 
have yielded crystalline structures with densities far exceeding that achievable using 
lithographic procedures [HO [3]. 

Much of the modern modeling of dissipation in systems of moving nano-particles is 
based on the work of George Gabriel Stokes, who formulated his famous drag law for 
the resistance of spherical particles moving in a viscous fluid at low Reynolds numbers. 
Reynolds number is defined as the ratio R = UL/u where U is the typical velocity, L 
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is the typical size and v is the kinematic viscosity. Reynolds number is physically the 
ratio between a typical inertial term and viscous force, and R -C 1 means the total 
dominance of viscous forces. It is commonly assumed that all processes in fluids at 
micro- and nano-scales are dominated by viscous forces and the Stokes approximation 
applies. The Stokes result states that a round particle moving through ambient fluid 
will experience resistance force that is proportional to the velocity of the particle [I]. 
Conversely, in the absence of inertia, the velocity of a particle will be proportional 
to the force applied to it since resistance force and applied force must balance. This 
law, that "force is proportional to velocity" is sometimes called Darcy's law 
Generalizing this concept to apply to an assembly of a large number of particles, we 
define the relative density of particles as the number of particles per unit volume, p. 
Given the energy of particles in given configuration E[p], the local velocity of particles 
under the force F = —V5E/5p is u = —pijFj. The coefficient of proportionality 
is called the mobility matrix p,^. For round particles, pij = pd~ij where d~ij is the 
Kronecker delta (equal to 1 if i = j and otherwise), so local velocity u is simply 
proportional to the local force F, u = pF. In most cases, the mobility is extended 
to depend on p itself. In particular here we allow for p to be an average value of the 
dynamical variable, which is given by a convolution p = K * p with some kernel K. 
The equation of motion for the density p is then given by the following conservation 
law whose general setting can be traced back to Debye and Hiickel [5] 

dp t-,5E 

— = -divpu = divp/iV— . (1) 

This equation neglects the diffusion of p which may occur if the particles also ex- 
perience Brownian motion. Note that the rate of change of the energy E[p] and an 
arbitrary functional of density F[p] is given by: 



dE f r ,f^SE\\ dF f . .SE 8Fj 



(2) 



which can be seen by direct application of equation (CQ). See also [HIE1E] for relevant 
discussions. 

Despite apparent mathematical differences, the work of Gilbert on dissipation in the 
Landau-Lifshitz equation for spin waves may be recast into the Darcy law form. 
Indeed, Gilbert [10] assumed a dissipation 'force' F[M] for a set of spins M to be 
proportional in the rate of change of M that we shall denote M: 

F$M\=VijMj. (3) 

Gilbert proposed that in the absence of any information on the matrix rjij it is prudent 
to assume simply rjij = f]Sij, which ultimately leads to the Landau-Lifshitz-Gilbert 
dissipation term [TU1 [TT] which we write as 

^ = 7 M x H+ ArM x M x H, (4) 

dt M 2 ' w 
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where H is the magnetic field, which can be expressed in terms of the external field 
H ex t and the energy of particle interactions i?[M] as H = H ex t + 5E/8M. (the varia- 
tional, or Frechet derivative) and 7, A are constants. We shall note that alternative 
phenomenological theories of dissipation exist (Bloch-Brombergen |13j . COT [14] ) 
that are not in a generalized Darcy's form. The Hamiltonian part of this equation 
(the first term on the right hand side) has the same dynamics as the rigid body, in 
terms of angular momentum. The important feature of this model is that it con- 
serves the magnitude of the magnetization vector |M|, while dissipating the total 
energy 7 J M ■ H d 3 cc. 

If we only take into account the dissipative (right hand side) terms in (J3J), then the 
rate of change in the energy i?[M] and an arbitrary functional of the magnetization 
M is given by 

dE f / 5E\ 2 , dF f / 5E\ / 6F\ , 

— = -/ Mx dx , — = - / M x • ( M x dx , (5) 

dt J V SMJ ' dt J V 5MJ \ dMJ ' v ; 

Thus, in essence, both the Stokes drag law and Gilbert dissipation are based on a 
single mathematical assumption: a 'dissipative force' acting on a physical quantity 
is proportional to the rate of change of that physical quantity. The Stokes drag law 
is valid for mobile particles that interact via collisions; Gilbert's work is developed 
for particles with the immobile center of mass that can orient themselves and change 
their magnetic properties. As far as we are aware, the work since Gilbert and Landau- 
Lifshitz on magnetic dissipation has remained largely phenomenological. A theory of 
magnetic dissipation for freely moving particles that is applicable to the nano-particle 
self-assembly still remains to be derived. The challenge is place the phenomenological 
arguments that have been shown historically to be practical for magnetic dissipation 
into a systematic theoretical framework. 

This paper derives equations for the dynamics of magnetized particles moving through 
a passive ambient fluid. Our derivation of the equations motion is not based on phe- 
nomenological arguments. Equations (T2J) and (0) possess intriguing physical and 
mathematical similarities, but a direct generalization to moving particles with ori- 
entation does not seem feasible. Instead, we proceed by introducing dissipation into 
the kinetic (Vlasov) equation describing the motion of particle density f{p,q,t) in 
the momentum p and coordinate q space. Only that approach can be physically 
justified. Our starting point is the formulation of the dissipative kinetic theory of 
gases consisting of arbitrarily shaped microscopic rigid bodies by accurately taking 
into account the linear and angular momentum of a microscopic particle. The con- 
tinuum equations of motion for the density and magnetization are then obtained by 
integrating out the momentum coordinate and truncating the expansion to the first 
order (a procedure known as 'moment truncation' in plasma physics). Remarkably, 
our general framework recovers both the continuum conservation law ([T]) for particles 
without orientation, and the Gilbert dissipation (J4]) for the oriented particles with the 
immobile center of mass. Moreover in some particular systems, we observe the for- 



5 



mations of domains of coherent orientation and we analytically predict the evolution 
of each domain. 

We note that an alternative form of dissipation in continuum media, caused by Brow- 
nian motion and expressed as diffusion terms in the Fokker-Planck equations is not 
suitable for our purposes, as applied to density, orientation or Vlasov form. The 
Fokker-Planck equation is a drift-diffusion extension of the Vlasov equation for ki- 
netic dynamics that is widely used as a powerful tool in understanding the physics of 
certain dissipative phenomena, especially in plasma physics. In this regard we men- 
tion some important results by Chavanis et al. (see [T3] and references therein), where 
a particular form of equation ([I]) with \x = const, has been obtained from a general- 
ized Fokker-Plank treatment in the simpler case of isotropic interactions. However, 
informed by recent results in nano-scale self-assembly, we assume that Brownian par- 
ticle diffusion may be neglected for our purposes. For other experiments with strong 
diffusive effects, one would seek an alternative theory. 



2 A dissipative kinetic equation for collisionless 

systems 

Recently, it was demonstrated [TB] that the Darcy's conservation law flTJ can be ob- 
tained from the introduction of dissipation in the kinetic Vlasov equations, describing 
the evolution of the density f(p,q,t) in the phase space: 

Inertial part Dissipation part 

where H[f] is the Vlasov Hamiltonian and 5H/5f is its variational (Frechet) deriva- 
tive, governing the single particle motion. The mobility (in phase space) fi[f] is a 
modeling choice and the energy functional E[f] is usually taken to be the potential 
part of the Hamiltonian H, describing just the self-interactions between the particles 
in the momentum space. If the inertial part is neglected, then by integrating this 
equation with respect to the momentum (p) we obtain precisely equation (TjQ), with 
p(q^) = / f(PiQit)dp an d fi\p] = J p[f(p,q,t)]dp. A particular form of equation (jSJ) 
with fi[f] = f was first suggested by Kandrup [17] in the context of galactic dynamics, 
and then analyzed later by Bloch et al. [HI [19] . The structure of this equation is 
consistent with the work of Kaufman [20J and Morrison [21J, where it is shown how 
dissipation can be introduced in Hamiltonian systems through a symmetric bracket, 
instead of the antisymmetric bracket, which is typical of Hamiltonian dynamics. This 
work was later extended by Grmela [22l [23| [23] . 
We list a series of relevant results that were proved in [16] . 
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Lemma 2.1 Equation (0|) preserves the single-particle solutions 

f = ^2w k {t)S(p-p k {t))S(q-q k (t)) , 

k 

(with tbk = 0) that are not allowed in other approaches in kinetic dissipation. 

Lemma 2.2 Equation (0|) preserves entropy S = J f log fdpdq and has infinitely 
many Casimirs: for an arbitrary function 4>{f), the integral J (p(f)dpdq is constant, 
independently of the choice of interaction energy. 

Lemma 2.3 If the inertial part is discarded, then truncating the moment equations 
and keeping only the terms dependent on p = J f(p,q,t)dp (zero-th moment terms) 
yields exactly Debye-Huckel equation (Qp. 

Let us now demonstrate how the dissipation given by ([6]) is connected to the Darcy's 
law. We neglect the inertial part of ([1]) and consider the rate of change for the 
energy E[f] and an arbitrary functional F[f] using Kandrup's choice for the mobility 

M = f- 

1 £--/{'.?}{'■£}*■ < 7 » 

This result can also be easily generalized to mobility functionals having the form 
fi[f] = f M[f], with M[f] > being a scalar functional of / [IE]. For general 
mobility functionals fi[f], the quadratic form given by energy dissipation functional 
(J7|) is not necessarily positive definite. 

Thus, just as equation ([I]), the dissipation in (jSJ) enforces the rate of change of an 
arbitrary functional of the solution to be linear in both 5E/5f and SF/Sf. In par- 
ticular, energy dissipation is given by a quadratic form of 5E/5f. This formulation 
in terms of energy dissipation is actually equivalent to saying 'force is proportional 
to velocity', but it has much more physical and mathematical power, which we are 
going to exploit in this paper. 

The mathematical equivalence of the dissipation laws (J2j [5]) and (JTj) is now clear: 

Assumption 2.4 The energy dissipation of an arbitrary functional F[f] is a quadratic 
form in the variational derivative of F[f] and the energy E[f] - see (TJj) below. 

There seems to be no direct path to consistently generalize the evolution equation 
for moving magnets to include both density and orientation, somehow combining (fTl) 
and (111) . The path we follow is to extend the kinetic equation ([6]) to account for 
orientation interaction and then compute the evolution equations using the moments 
approach. 
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3 Historic digression and mathematical background 

3.1 Selective decay hypothesis 

The inspiration for this work originated in the dissipative bracket structure intro- 
duced in Bloch et al. pj5] , which in turn was motivated in part by the double bracket 
introduced by Vallis, Carnevale, and Young [28] for incompressible fluid flows. The 
dominant idea for [28] was the selective decay hypothesis, which arose in turbulence 
research [29J and is consistent with the preservation of coadjoint orbits. According 
to the selective decay hypothesis, energy in strongly nonequilibrium statistical sys- 
tems tends to decay much faster than certain other ideally conserved properties. In 
particular, energy decays much faster in such systems than those "kinematic" or "ge- 
ometric" properties that would have been preserved in the ideal nondissipative limit 
independently of the choice of the Hamiltonian. Examples are Casimir functions for 
the Lie-Poisson formulations of various ideal fluid models [30] . 

The selective decay hypothesis was inspired by a famous example; namely, that en- 
strophy decays much more slowly than kinetic energy in 2D incompressible fluid 
turbulence. Kraichnan [31] showed that the decay of kinetic energy under the preser- 
vation of enstrophy causes dramatic effects in 2D turbulence. Namely, it causes the 
well known "inverse cascade" of kinetic energy to larger scales, rather than the "for- 
ward cascade" of energy to smaller scales that is observed in 3D turbulence! In 2D 
ideal incompressible fluid flow the enstrophy (the L 2 norm of the vorticity) is pre- 
served on coadjoint orbits. That is, enstrophy is a Casimir of the Lie-Poisson bracket 
in the Hamiltonian formulation of the 2D Euler fluid equations. Vallis et al. [28] 
chose a form of dissipation that was expressible as a double Lie-Poisson bracket. This 
choice of dissipation preserved the enstrophy and thereby enforced the selective decay 
hypothesis for all 2D incompressible fluid solutions, laminar as well as turbulent. 
Once its dramatic effects were recognized in 2D turbulence, selective decay was posited 
as a governing mechanism in other systems, particularly in statistical behavior of fluid 
systems with high variability. For example, the slow decay of magnetic helicity was 
popularly invoked as a possible means of obtaining magnetically confined plasmas |32j . 
Likewise, in geophysical fluid flows, the slow decay of potential vorticity (PV) relative 
to kinetic energy strongly influences the dynamics of weather and climate patterns 
much as in the inverse cascade tendency in 2D turbulence. The use of selective decay 
ideas for PV thinking in meteorology and atmospheric science has become standard 
practice since the fundamental work in [33l IM] . 

3.2 Mathematical framework for geometric dissipation 

As explained in [7], dissipation of energy E may naturally summon an appropriate 
metric tensor. In previous work, Holm and Putkaradze [351 EHj showed that for any 
two functionals F[k], G[k] of a geometric quantity k a distance between them may be 
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Here ( • , ■ ) denotes the L 2 pairing of vector fields (X) with their dual one-form densi- 
ties (X*), sharp ( ■ )" raises the vector index from covariant to contravariant and //[k] 
is the mobility functional. The mobility //[k] is assumed to satisfy the requirements 
for (jSJ) to be positive definite and symmetric, as discussed in [39]. The diamond op- 
eration (o) in equation (jSj) is the dual of the Lie algebra action, defined as follows. 
Let a vector field £ act on a vector space V by Lie derivation, so that the Lie algebra 
action of £ on any element k G V is given by the Lie derivative, 

The operation dual to the Lie derivative is denoted by and defined in terms of the 
L 2 pairings between X and X* and between V and V* as 

Given the metric (JSJ) and a dissipated energy functional E[k], the time evolution of 
arbitrary functional F[k] is given by [35j [39] as 

f = flF , S»M := - 9 .(F,B) = - ((„[«] o f ), (« o f )")^ , (10) 

which specifies the dynamics of any functional F[k], given the the energy dependence 
E[k\. The bracket {{ F , E }} is shown to satisfy the Leibnitz product-rule property 
for a suitable class of mobility functionals in [351 [39]. Eq. ffTUj) and positivity of 
g K (E, E) imply that the energy E decays in time until it eventually reaches a critical 
point, 5E/5k = 0. 

Remark 3.1 For densities (dual to functions in the L 2 pairing), the Lie derivative is 
the divergence and its dual operation is (minus) the gradient. Thus, for densities the 
symbol diamond (o) is replaced by gradient ( V) in the metric defined in Eq. l[TU\) . 

4 Kinetic equation for anisotropic particles and 
Landau-Lifshitz equations 

We now turn to the problem of a system composed of immobile particles of arbitrary 
shape. The orientation m of particles now appears as an additional degree of freedom, 
but the particle momentum p is going to be neglected in this section since centers 
of mass for the particles are assumed to be immobile. In our approach we consider 



(9) 
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orientation as given by the angular momentum of each particle. Thus, instead of 
considering unit vectors, we consider angular momentum vectors and each particle is 
treated as a rigid body. Let us consider the non-dissipative case first. The evolution 
equation for the distribution function <p(x,m,t) in this approach is written in the 
conservation form in the U = (U x , U m ) on the (x, m) space 



= -div (a!i m) O U) = - Va, ■ U a 



d 
dm 



As we mentioned, in this section we set U x = since the centers of mass for the 
particles are assumed to be stationary. Since the angular momentum for a rigid body 
evolves according to the Hamiltonian motion m = m x V m /i = U mi then we can 
write our kinetic equation as 



dip 



5H 

5<p 







where we have used the definition of the rigid body bracket [tp, h] := m ■ V m <£> x 
V m h and the global Hamiltonian is written as H[<p] = JJ (p(x, m) h(x, m) d 3 x d 3 m. 
Now that we have obtained the conservative dynamics in the bracket form, we are 
motivated to insert the dissipation as the double bracket form of §5§ as follows: 



dip 
~di 



5H~ 








, 8E~ 




^' 5p 








^Jp~_ 





Inertia 



Dissipation 



In this context, it is reasonable to take the moments with respect to 
lar we want the equation for the magnetization density M(£C, i) := 
H = E. We now introduce the assumption that the single particle 
in m (it is easy to see that this assumption recovers, for example, 
the potential energy of a magnetic moment in a magnetic field). 
Again, let us consider only the dissipation terms in (jlip . i.e., the 
hand side and use Kandrup's mobility fj,[ip] = ip. Then, we get 
energy E[tp] as well as an evolution of an arbitrary functional F[i 
form analogous to 



(11) 



m and in particu- 
jrrup d 3 m, when 
energy h is linear 
the expression of 

term on the right 
the dissipation of 
] in the quadratic 



dE r 


5E~ 


~dT~ ~J 


^' 5(p 



dx 



dF 
~dt 



5E 

S(p 



6E 

5<p 



dx. 



(12) 



The extension to more general mobility functional is also possible, but the quadratic 
form describing the dissipation will not be symmetric except for some particular 
choices for the mobility. 

Now, using the method of moments we compute the equation for the magnetization 
density M as 

dM 5H _ SE 

M x h M x fi M x 



dt 



5M 



5M 
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recovering exactly the Landau-Lifshitz-Gilbert equation (LLG) for magnetization dy- 
namics (J4]) for /z[M] = M, which occurs for Kandrup's original choice of mobility, 

M = /■ 

It is interesting to note that similar ideas for the extensions of Gilbert dissipation 
were recently suggested phenomeno logically in [36] . 



Dissipative kinetic equation for moving particles 
with orientation 



The kinetic approach for immobile oriented particles can be generalized to moving 
oriented particles. We are using the energy dissipation approach by generalizing the 
formula ([7J to more general Poisson bracket that involves both density and orienta- 
tion. Interestingly enough, this Poisson bracket has been derived 25 years ago in the 
[371 EH] > but we believe that we are the first to apply this bracket to dissipation in 
this setting. This modified Poisson bracket {■ , -}i is defined as 

/' v}i = r*7. 

that is, it is a sum of a canonical Poisson bracket {■ , ■} and the Rigid body bracket 
[• , ■]. For g belonging to the dual of a general Lie Algebra g (not necessarily so(3)*), 
we can compute elements SE/5g G Q and SF/Sg e q. We thus define a generalized 
rigid body bracket as 

~5E 5F~ 



+ 



6E 



(13) 
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6g ' 5g 



LA 



(14) 



Here, we have used [, ]la ~ the Lie algebra bracket in q. For this choice of the new 
bracket the quadratic form describing the dissipation (J7J generalizes naturally to the 
anisotropic follows: 



dF 
~dt 



5E 



(15) 



Taking the moments of the anisotropic kinetic Vlasov equation following from ffTSI) 
leads to a new kind of continuum dynamics naturally generalizing equation ((TJ) to the 
case of anisotropic interaction. To simplify the formulas, we neglect the inertia of the 
particles. As mentioned in the introduction, this is a common assumption for friction 
dominated systems. This assumption means neglecting the term {/, SH/5f} in the 
left hand side of Thus we are left with the following equation which will form 
the core of further analysis, 



df 
dt 



f,<M, 



5E 



Sf 1 



(16) 



The remainder of this section follows the reasonong in [39J , where the reader may find 
many more technical details and exposition of some related problems. 
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5.1 Landau-Lifshitz through the Smoluchowski approach 



When dealing with kinetic theories, the kinetic moments are a well established tool 
for deriving simplified continuum descriptions. However, in the case of anisotropic 
interactions, the kinetic moments A n = J p n f(q,p, g) dp depend not only on the 
position coordinate, but also on the single-particle Lie algebra coordinate g. Thus, 
the moment equations differ from those found in the isotropic case, since they must 
include the terms which take into account the anisotropic behavior. 
Moment dynamics has a Hamiltonian structure which is justified by the theorem 
that taking moments is a Poisson map. In fact, their dynamics is governed by the 
Kupershmidt-Manin bracket [26j EZ] • This process can be transferred to the double 
bracket formulation as follows. As in the Kupershmidt-Manin approach, the moments 
are dual to the variables /3 n (q, g), which are introduced by expanding the Hamiltonian 
function h(q,p,g) as h(q,p,g) = p n fl n (q,g). The Lie algebra action is given by 

Pn-f = {f,p n Pn} l 

The dual action is defined by a star (*) operator which is the symplectic equivalent 
of the diamond (o) operation, as follows, 

(/*„*, (3 n ) := (/, n k) = (f*k,p n (3 n (q,g)) = /f p n {f,k} l( ip, (3 n 



Thus, the star operator is defined explicitly for k = p m a m as 

dA m+n da m 



f ~k n k — ad Am+n-i 



9- 



dg dg 



where the ad* operator is the Kupershmidt-Manin Hamiltonian operator, see [261 12?] . 
We introduce the dissipative bracket by 



F = {{F,E}} 



(17) 



At the simplest level we truncate this sum, by taking n — 0. For higher order 
approximations, giving highly complex formulas, see [39]. By using this evolution 
equation for an arbitrary functional F, the rate of change for zero-th moment A is 
found to be 



dt 



fj- 



9A <9tq 
dg ' dg 



where one defines 70 := fi[f] *o 5E/5f, which is explicitly given by 

~dp d/3 ' 



To 



dg ' dg 
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Upon introducing the notation {•, •} for the Lie-Poisson bracket on the Lie algebra g 

{k, h} :-- 



g, 



dk dh 
dg : dg 

we may write the equation for ip = A in the double bracket form 



dip 



5E 



(18) 



This is the purely dissipative moment equation corresponding to the Vlasov equation 
for the oriented particles ffTBT) . It is straightforward to see that the inclusion of the 
inertial part such as in equation yields a moment equation of exactly the same 
form as equation (fTTT) 



dip 
~dt 



SH_ 

5<p 



5E 

5ip 



(19) 



We have now shown how to derive equation from the kinetic approach by taking 
the zero-th moment. Next, one could continue by integrating over the microscopic 
magnetization m and thereby obtain the Landau-Lisfshitz-Gilbert equation (jlj). 
In the particular case g = m £ so(3) that is of interest to us, the Lie-Poisson bracket 
reduces to the rigid body bracket {ip, h} := m-V m ip x V m h. The global Hamiltonian 
is written as H[ip] = ffip(q,m)h(q,m)dqd 3 m and an analogous expression holds 
for the dissipation energy E[ip\. The generalization to higher spatial dimensions is 
straightforward, since there are no derivatives in spatial coordinates in the equation 

We next take moments with respect to m. The particular physical case of magne- 
tization density M(cc,i) := Jmipd 3 m requires H = E. We assume that the single 
particle energy h is linear in m since this assumption recovers the expression of the 
potential energy of a magnetic moment in a magnetic field. We consider only the 
dissipation terms in (fTTj) . i.e., we drop the second term term on the left hand side of 
that equation. Then, we obtain the rate of change of the energy E[ip], as well as of 
an arbitrary functional F[ip\ in the quadratic form analogous to ([7j): 



dE 
~dt 



dx ' H 



' U] 5E 



f ~ 



dx . (20) 



The method of moments gives the equation for the magnetization density M 



dM 
~dt 



5H ^ 6E 
M x + M x /Lt[M] x 



5M 



5M ' 



which recovers exactly the Landau-Lifshitz-Gilbert equation (LLG) for magnetization 
dynamics (j4j) for /x[M] = M, which occurs for Kandrup's original choice of mobility 

M = /■ 
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5.2 Dissipative moment dynamics: an alternative treatment 

In this section, we extend the Kupershidt-Manin approach as in [37J [38] to generate 
the dynamics of moments with respect to both momentum p and charge g. The main 
technical complication is that the Lie algebras of physical interest (such as so (3)) 
are not one-dimensional and in general are not Abelian. Thus, in the general case 
one needs to use a multi- index notation as in [251 126]. We introduce multi- indices 
a := (<7i, 02, . . . , cr N ), with cr, > 0, and define g a := g* 1 . . . g^f , where N = dim(g). 
Then, the moments are expressed as 

A na (q) := J p n g a f(q,p,g)dpdg. 

This multi-dimensional treatment encumbers the calculations with technical complli- 
cations. For the purposes of this section, we are primarily interested in the equations 
for p and G, so we restrict our considerations only to the moments of the form 

An,v = J P n 9uf(q,P,g) dpdg v = 0, 1, ... ,7V. 

Here we define go = 1 and g a = (g, e a ) where e a is a basis of the Lie algebra and 
(gbG b , e a ) = g a £ K represents the result of the pairing ( ■ , ■ ) between an element 
of the Lie algebra basis and an element of the dual Lie algebra. We write the single 
particle Hamiltonian as h = 5H/5f = p n g u 5H / 5A ritU =: p n g u (3 n '(q), which means 
that we employ the following assumption. 

Assumption 5.1 The single-particle Hamiltonian h = 5H/5f is linear in g and can 
be expressed as 

h{q,p,g)=p n i) n {q)+p n {g^ n {q)) , 

where ip n (q) £ R is a real scalar function andi/j n (q) £ Rxq is a real Lie- algebra-valued 
function. 



Dual Lie algebra action. The action of (3^ on / is defined as 

P:-f={p n 9vP:, f} x (no sum). 

The dual of this action is denoted by (*„,„). It may be computed analogously to Q 
and found to be 

/ *n,u k = J J p n g u {/, k\ <\pdg 

.* a i , / / \ dA m+n d(g a a a m )]\ 
g u g a ad a a n A m+n _idg + J g„(g, — , — — ) dg 

A* f A A , f I \ 9A m+n d(g a OL%) 1 \ 

ad„, CT / g v g a A m+n _i d^f + / g„(g, — ^ , ~ ) &9 ■ 



14 



Here, k = p m g a ac^q) and we have used the definition of the moment 



A n (g,g) = J p n f(q,p,g)dp- 

Evolution equation. Having characterized the dual Lie algebra action, we may 
now write the evolution equation for an arbitrary functional F in terms of the dissi- 
pative bracket as follows: 



F = {{F,E}} 



(21) 



where the pairing ((■,■)) is given by integration over the spatial coordinate q. Now 
we fix m = 0, n = 1. The equation for the evolution of F — A 0t \ := f g\A dgdp is 
found from (ED) to be 



dA , 



x 



Of 



^ f 7i,„ J 9v 9\ Aq dg^j 



9x { 9, 



dAt d{g al l a ) 



dg dg 



dA d{g al l a ) 



dg dg 



where we have defined the analogues of Darcy's velocities: 



7o,^ : 
and 



dg ' dg 



dg 



9u \ 9, 



dg 



d^o d(g a f3%) 



(22) 



71,1/ := n\f\ 



5E 



Sf dq 



gug a ^dg+ I g u (g 



dg dg 
dg ' dg 



dg 



dg. 



with Jlk '■= J p k f dp. Here we have assumed that the energy functional E depends 
only on Aq^\ (recall that (3^ := 5E/5A n± \), so that we may fix k = in the second 
line. At this point, we further simplify the treatment by discarding all terms in 7i j0 , 
that is we truncate the summations in equation ( J22l) to consider only terms in 70,0, 
7o, a and 71,0. With this simplification, equation (1221 becomes 



OA, 



d 



dt dq 
and the expression for 7i j0 is 



7i,o / 9x A dg + / g x ( g, 



dA d(g alo l) 



dg dg 



dg 



(23) 



7i,o := ml *i,o jj = -q^- I 9a Vodg. 
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Let us now define the dynamic quantities 

P = //d S dp, G = fgfigip. 

and corresponding mobilities 

Hp = J n[f]dgdp, fiQ = J gp[f}dgdp. 

In terms of these quantities, we may write the following 
Theorem 5.2 The moment equations for p and G are given by 
dp d ( ( d 6E I d 5E\\ 

m = Y q { p [^d- q Jp- + \^d- q lG/) I (24) 

and 



dG d 8 8E I dSE\\\ 



Remark 5.3 Similar equations (under the name of Geometric Order Parameter equa- 
tions) were derived via a general geometric considerations in [35^ . The equations 
presented here reduce to those for the case of commutative Lie algebras. 

Remark 5.4 Although in this section the underlying space where particles undergo 
their motion is assumed to be one dimensional, the extension to higher spatial dimen- 
sions is trivial: since the Lie derivative of a density is given by the divergence term 
in spatial coordinate q ; we need to substitute div and grad instead of one- dimensional 
derivatives with respect to q. See ft2hj27\ ) immediately below for an explicitly computed 
particular case. 



6 Self-organization dynamics of moving magnets 

We are now ready to write out explicit equations for the motion of magnets in d = 1,2 
or 3 dimensions. In this case G — M(a?) Gl 3 - average magnetization, ad v w = v x w 
and ad* w = — v x w, and the Lie algebra pairing is represented by the dot product of 
vectors in IR 3 . Therefore, the equations for (p, M) are obtained directly from (I24|25l) 
and are written as follows: 



(26) 
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and 



dM 



dt 



div M <8> ( A^pV— + n m ■ V 



5M 



5E 



) 



+ M x n M x 



5M 



5E 



(27) 



This is the generalization of the Landau-Lifshitz-Gilbert and Debye-Hiickel equations 
we have sought. The main consequence of equations f|2TI) is that they allow single 
particle solutions of the form 



These solutions have an important physical meaning, since they represent particles 
that aggregate and align. This phenomenon is of fundamental importance in the the- 
ory of anisotropic self-assembly. There are special cases where these solutions emerge 
spontaneously from any confined initial distribution. These behavior depends on the 
particular choice of mobilities and energy functionals. A typical result of simulations 
is presented on Fig. [TJ Starting with a random distribution in one dimension, we see 
the formation of sharp peaks in magnetization M (left figure) and density p (color 
code on the left and plot on the right), corresponding to the (^-functions. 

7 Moving dissipative oriented curves: new vari- 
ables in geometric dynamics 

The structure of singular solutions allows an interesting geometric generalization, 
namely a possibility for singular solutions to be concentrated on a one- dimensional 
manifold (a curve) in the two- or three- dimensional space. This will give an interest- 
ing addition to the theory of exact geometric rods pEO] . as our equations describe the 
dissipative motion of a curve. Formally, we extend equations ff28l) and assume that 
the solution is concentrated on a curve in space r = Q(s, t), where s is the arclength. 
Then, we are interested in the solutions of the form 



p(x,t)=w p (t) 5(x-Q(t)) 
M(x,t) =w M (t) 5{x-Q{t)) 



(28) 



where w p , u>m and Q undergo the following dynamics 




p(x,t) = w p (s,t) 6(x - Q(s,t)) 
M(x,t) = w M (s,t) S(x - Q(s,t)) 



(29) 
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t=0 t=2 




Figure 1: Left: A typical example of orienton formation in a d = 1 dimensional 
simulation. The color-code on the cylinder denotes local averaged density ~p = H * p 
for a Helmholtz kernel H = e~' x ': black is maximum density while white is p = 0. 
Purple lines denote the three-dimensional vector M = i?*M. The formation of sharp 
peaks in averaged quantities corresponds to the formation of ^-functions in ( 1261 1271) 
according to (l28il . Averaged quantities were chosen to avoid the necessity to represent 
(^-functions. Right: The corresponding waterfall plot for evolution of averaged density 
J> = H * p. Horizontal coordinate is space. Sharp peaks correspond to the formation 
of (5-function singularities in the density variable p. 
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Multiplying each equation fl26| 12?]) by a test function and integrating produces the 
following equations of motion: 

J^p(M) = o (30) 

— w M (s,t) = w M x n M (Q(s,t)) x _(Q( s ,t)) (31) 

J^Q(M) = ti p (s,t)V^(Q(s,t)") +AlM ( s ,t)V^(Q(s,t)) (32) 
In the simplest case of a binary energy functional, 

E = ^J p(x)p(x')G P (x - x') + M(x)M(x')G M (x - x') d 3 x d V , (33) 

where Gm is assumed to be a scalar function, equation (132!) assumes an especially 
simple and elegant shape - see also [JT]: 

-Q( S ,t) = J ii p {s,t)p{s',t)VG p (Q{s,t),Q{s',t)) 

+ /x M (s, t) • M( S ', t) VG M (q(s, t), Q(s', t)) ds' . (34) 



7.1 Numerical solution of curve evolution 

To illustrate the concepts developed above, we perform a simulation of the spatio- 
temporal evolution for two oriented curves, having the simplest energy dependence 
given by (133]) . In the future, we shall consider more sophisticated models of energy 
using the technique developed in ( l47j) below. In computer simulations presented 
on Fig. El we assumed centrally symmetric Lennard- Jones potential for the density 
interaction and the Helmholtzian Gm(x) = exp(— |x|) for the orientation part of 
energy. Initially, the curves are far away from each other and the M values on the 
curves are unrelated. Eventually, the curves collapse to two parallel lines and the M 
vectors on both curves align. 



8 Energy dependent on interactions of orientations 

Energy (133]) provides insight, but it is too simple to describe many interesting physical 
phenomena, for example, dynamics of assembly of biologically relevant particles into 
curves like DNA. The energy of these particles depends on the position x and orien- 
tation A(x) of a frame at a point x, and the relative distances and orientations with 
respect to all other points x'. We are thus compelled to deal with more interesting 
energies reflecting this physical fact, and this will lead us to some interesting mathe- 
matical consequences. The technical question here is how to compute 5E/5M., if the 




Figure 2: An example of two oriented curves (red and green) attracting each other 
and unwinding at the same time. The blue vectors illustrate the vector M £ so(3) at 
each point on the curve. Time scale is arbitrary. 
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energy is given in terms of A. Of particular interest to us is the so-called anisotropic 
Lennard- Jones potential that has been developed in [12] and has been subsequently 
applied to produce self-organization of anisotropic particles into the shape of a double 
helix [13]. The idea of the method is the following. 

Consider two ellipsoids with the geometric centers at x and x' and orientation with 
respect to the fixed frame A(x) and A(x'), respectively. The method computes the 
distance d(x, x', A(x), A(x')) between ellipsoidal surfaces along the direction of x — x'. 
This is close to the minimal distance between the ellipsoids, but it is easier to compute. 
In principle, to make the depth of the potential well different for attracting and 
repulsive potential, one can consider two distances d — d\ and d = d 2 , corresponding 
to the attractive and repulsive parts of the ellipsoids. The distances d± and d 2 are 
defined as di = <To + |x— x'|+/i(A(x, Ax')), where o"o is a given constant and /, are some 
(rather complex) functions. One then takes the potential to be the Lennard- Jones 
potential of the modified distances: 

t/{rf(x,x / ,A(x),A(x'))} = 4e 

An important feature of this potential that has not been noticed in previous works is 
its invariance with respect to 50(3) rotations about each point on the curve. Thus, 
the energy functional may be written as 

E = \J p(x)p(x')t/(x,x', A(x), A(x'))dxdx', (36) 

where U is some scalar function of its arguments. This energy is still binary, i.e. 
it is a direct sum of pairwise interaction of all particles, which is known to be an 
rather crude approximation for real systems. However, we believe it is important to 
first understand this 50 (3)-invariant energy. Later, the equations of motion for more 
general energies may be obtained by direct generalization. Invariance with respect 
to the rotation group requires that the energy of particle interaction assumes the 
following form: 

E = \j p(x)p(x')f/(|x-x / U(x,x'))dxdx', (37) 
where we have defined 

£(x, x') = A(x)- 1 A(x / ) e 50(3) . (38) 

Remark 8.1 Note that we have assumed that the energy does not depend directly on 
a Lie algebra element M = A~ X VA. That means, in particular, that there is no elastic 
energy caused by particle motion. The introduction of elastic energy complicates the 
equations and will be considered in future work. 



(70 

di 



12 



d 2 - 



(35) 
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Now, we need to compute 5E/SA. We proceed as follows. Define £(x) = A 1 (x)5A(x) e 
so(3). Then, notice that 

= -A- 1 (x)M(x)A- 1 (x)A(x / ) + A^ 1 (x)M(x / ) 

= -E(x)£(x,x') + A- 1 (x)A(x / )A- 1 (x / )M(x / ) 

= -£(x)£(x,x')+£(x,x')X(x'), (39) 

since £ = A" 1 (x)A(x / ) e 50(3). Thus, 

SE= |p(x)p(x')tr^^(x,x')) {-S(x)^(x,x')+e(x,x')S(x')}jdxdx' 

= | p(x)p(x')tr (- ^(x,x')e(x,x') T ^S(x)j dxdx' 

+ | p(x)p(x')tr ^e(x,x') T ^(x,x')) T S(x')j dxdx'. (40) 

In the last term, x and x' were exchanged to make the variation S(x') a function of 
x. When making this exchange, we must remember that 

£(x,x') = A(x)" 1 A(x') = T^x'.x) = £ T (x',x) , 

and similarly, 

Consequently, we find, 

^ = / P(x)p(x')tr (- ^(x,x')e(x,x') T ) T S(x) 

+ | p(x)p(x')tr j^(x,x') (^(x,x'))^ S(x)J dxdx'. 

Thus, collecting terms leads to 

5E = J p(x)p(x')^-^(x,x')e(x,x') T + e(x,x') ^(x,x')) T ,S(x)^)dxdx'. 

If we now turn from the variations in £ to variations in A, we find A in front of the 
integral, and the final formula for the variation of energy is thus 

^ = A(x)p(x) |p(x'){-^(x,x'r(x,x')+e(x,x')(^(x,x')) T }dx'. (41) 
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We have computed the variations with respect to A. We are not finished yet, as the 
equations of motion of a string are formulated on the Lie algebra, and not the tangent 
manifold to Lie group. Thus, we need to connect 5E/5A to <5£/5M0 To do that, we 
use the definition of the left-invariant tangent space at the identity, 

<9A 

M^A" 1 — eso(3), i = l,2,...,d, (42) 

where x l is the i-th coordinate in space. The Lie algebra element is then a 'vector' 
of d antisymmetric matrices, where d is the dimension of the space. We have, by the 
usual change of variables (summation over repeated indices is assumed), 

«-(£■">-(£•*>■ (43) 

where SE/SM e so(3)*. Denote again £(x) = A -1 (x)5A(x) G so(3). Then, 

Q 

5Mi = -EM, + A^—SA . 

ox 1 

Thus, ( 1131) becomes 



/ \M ^ / 



/ 5E 



\5Mi 



M, : , E 



9 5£ 



(44) 



Remembering that E = A 1 SA and setting A using the pairing properties 
< a, Ab >=< A T a, b > leads to the following expression (as (A~ X ) T = A): 

5E . { 6E d 8E\ 

— = A ad~ — =— - . (45) 

This can be written in a familiar way in terms of the vector product in M 3 as 

5E . ( -\ jr SE d 5E \ 

— = -A M ? - x h — , (46) 

which is an extension of the formula (6.18d) of [40] . 

Combining equations (146]) and (1411) yields an auxiliary equation for 5E/5~M. that is 
valid for any value of x: 

( Mi (x)x— (x) + ^ — (x)^ =-p(x) /'p(x , )(-?(x,x / )e(x,x') T 

+ ^(x,x')(^(x,x')) T }dx'. (47) 



1 Here, the hat denotes the usual isomorphism between vectors in M 3 and elements of Lie Algebra 
so(3), so M G so(3) and M € M 3 . 
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where, again, £(x, x') := A~ 1 (x)A(x / ) G SO(3). The left hand side of this equation 
is the so(3)*-valued covariant divergence of the so(3)* vector 5E/5~M.i, which is also 
a 3-by-3 antisymmetric matrix for each i = l,2,...,d. Likewise, the expression in 
curly brackets is an antisymmetric 3-by-3 matrix, as one may verify by taking the 
transpose. Thus equation (147j) is consistent. In fact, equation (1471) is reminiscent 
of the Yang-Mills Gauss Law for an so(3)*- valued gauge charge. A similar equation 
appears in the gauge-theoretical formulation of the dynamical equations for a spin- 
glass [H] . The gauge freedom remaining in equation (T4"T|) is removed by noticing that 
Mj defined in ( T4"2"j) is a pure gauge field. Consequently, the connection associated 
with the so(3)-valued covariant divergence in (T4T1) must have no curvature. Hence, 
for any i, j, we find the following consistency conditions 



dMi dM 



3 



Mi , M,- 



M;M ; - M i M ! , (48) 



dxi dx l 
which close the system ( l4Tj) . 

Remark 8.2 It is interesting to note the similarity between the variable £ = A~ 1 (x)A(x / ) 
and its time analogue A(t)~ l A(t+ A) , which is commonly used in discrete-time Moser- 
Veselov integrators JJ3j/ and their extensions U^l - Further study of the geometric 
nature of ( fj^l ) appears intriguing and will be undertaken in the future. 

Let us see how (T4T|) reduces on the singular solutions (1291) . Assuming that the poten- 
tial U is such that dU /d£ remains finite on the singular solutions, we substitute (1291) 
into (JH]) to find 

8E 

— (x) = A(x)^(s, t)5(x - Q(s, t)) 

J wptf,t){ ~ ^Mt T M+ZM (^( s ' s ')) T }dx', (49) 
where we have denoted 

r)JJ r)TT 

Z(s,s'):=t(Q(s,t),Q(s',t)) and ^±( s , s') := ^(Q(s,t), Q(s> ,t)) . 

Upon multiplying each side of the equation (14?]) by an arbitrary function 0(x), in- 
tegrating over the whole space and equating the coefficients of 0(Q(s, £)), we obtain 
the following relation: 



. , 5E 5E 

-m(,)x— (,) + - — =- Wp[ s 



as,s')(^(s,s')\ T }ds'. (50) 



This is the equation for the computation of the variations SE/SM. on the curve that 
we have sought. The analogue for curves of the consistency condition (1481) then closes 
the system. Future work will explore its dynamics and numerical solutions. 
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9 Conclusions 

In this paper, we explained the derivation of continuum equations for dissipation using 
the moment approach in kinetic (Vlasov) equations. We showed how to introduce 
dissipation in Vlasov equations based on a geometric generalization of the Darcy's 
law of motion (force being proportional to velocity). We derived the continuum 
equations of motions in general settings and then formulated those equations for the 
particular case of self-organization of magnetic particles. We showed the existence 
and spontaneous emergence of generalized solutions concentrated on delta-functions. 
We also showed that the equations possess singular solutions concentrated on one- 
dimensional manifolds (curves) and derived equations of motion for these curves. 
Of particular interest for future work are two questions. The first question is, Can so- 
lutions defining singular curves appear from random initial conditions in two or three 
dimensions? The positive answer to this question would demonstrate the possibility 
of controlling magnetic self-assembly of dispersed particles into "magnetic strings" , 
which would be very interesting technologically. 

The second question is, Can we consistently incorporate both inertia and dissipative 
terms in the theory of exact geometric rods? Hopefully, this can be accomplished by 
starting with the dissipative Vlasov equation containing the inertia terms, repeating 
the procedure of [37J EE] keeping the momentum term and closing the system using 
either truncation or cold plasma approximation. An interesting task would be to 
compare the resulting dissipative equations with those of the exact geometric rods in 
the absence of dissipation [40] . 
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